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This  paper  deals  with  the  extension  of  a crack  growth  model  to  high 
temperature  complex  loading  and  application  to  turbine  disc.  The 
proposed  model  is  based  on  an  extensive  experimental  study  performed 
by  SNECMA  Moteurs  on  N18  from  450°C  to  650°C,  in  isothermal  and 
also  in  non  isothermal  condition,  which  comprises  fatigue  with  or 
without  hold  times  and  special  sequence  tests  representative  to  the  disc 
in  service.  The  crack  growth  model  is  built  up  in  the  framework  of 
classical  linear  elastic  fracture  mechanics.  Time  effects  at  high 
temperature  are  traduced  by  creep-fatigue  and  oxidation-fatigue 
interactions.  The  proposed  writing  in  non  isothermal  condition  is  very 
attractive  for  easy  model  identification  on  a large  temperature  domain. 
Model  predictions  are  shown  for  a large  set  of  experimental  data 
including  complex  loading  in  non  isothermal  condition. 


1 INTRODUCTION 

This  paper  presents  a phenomenological  model  able  to  simulate  the  crack  growth  in  a 
metallic  material  at  high  temperature  and  under  complex  cyclic  loading  condition.  This  one 
uses  the  actual  stress  intensity  factor  as  the  loading  variable,  but  also  several  material 
parameters  are  introduced  in  order  to  describe  the  history  of  the  various  processes  that 
operate  near  the  crack  tip  during  for  instance  an  hold  period  at  high  temperature:  the 
mechanical  stress  redistribution  due  to  visco-plasticity,  the  creep  damage  process  under 
tensile  condition,  the  environmental  driven  brittle  oxide  at  crap  tip  that  increases  the  fatigue 
crack  propagation  and  the  interaction  of  overloads  with  both  fatigue  and  creep  crack 
growth. 

After  reminding  and  discussing  the  main  features  of  the  model,  originally  proposed  by 
Prigent  [1]  and  Kruch  et  al  [2],  we  propose  in  a first  section  a model  extension  to  account 
for  non  isothermal  complex  loading  together  with  the  functions  describing  the  fatigue, 
creep  and  oxidation  damages  and  the  different  interactions.  We  give  then  in  the  second 
section  the  results  of  the  application  made  with  N18  aeronautical  superalloy,  the  model 
predictions  of  many  isothermal  tests  performed  at  different  temperatures  and  also  of  non 
isothermal  tests  more  representative  of  real  loading  subjected  to  the  component  in  service. 


Paper  presented  at  the  RTO  A VT  Symposium  on  “Ageing  Mechanisms  and  Control: 

Part  B - Monitoring  and  Management  of  Gas  Turbine  Fleets  for  Extended  Life  and  Reduced  Costs  ”, 
held  in  Manchester,  UK,  8-11  October  2001,  and  published  in  RTO-MP-079(1). 
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2 CRACK  GROWTH  MODEL 

2.1  Assumptions  of  the  creep-fatigue-oxidation  interaction  under  complex  loading 

The  proposed  model  is  built  up  in  the  framework  of  classical  linear  elastic  fracture 
mechanics.  The  main  points  the  model  must  take  into  account  are:  the  fatigue  crack 
propagation,  the  creep  crack  propagation,  the  stress  intensity  factor  relaxation  during  creep, 
the  creep-fatigue  interaction  and  the  environment  effects.  The  various  state  variables  that 
store  the  post-loading  events  are  as  follows: 

- A threshold  for  fatigue  crack  growth,  similar  to  the  opening  stress  intensity  factor. 

Its  evolution  is  influenced  both,  by  fatigue  overloads  (and  under-loads)  as  proposed  by 
Baudin  and  Robert  [3],  and  by  the  hold  time  or  low  loading  rate  during  which  a 
relaxation  can  occur.  Built  up  in  the  basis  of  Onera’s  crack  propagation  model  proposed 
to  account  for  complex  loading  in  fatigue,  it  can  take  into  account  the  overload  effects 
applied  during  a cycle  through  the  plastic  zone  evolution  in  front  of  the  crack  and  the 
value  of  the  fatigue  threshold  K{ ; that  will  not  be  developed  in  this  paper. 

- This  threshold  is  also  introduced  to  describe  the  creep  crack  growth  Kcs  in  a time 
dependent  expression. 

- The  size  of  the  oxidised  zone  at  the  crack  tip  / , from  which  is  deduced  an  effective 
local  toughness  of  the  material  Kc , that  operates  in  the  fatigue  crack  growth  model 
through  a Forman’s  expression. 

To  account  for  the  creep- fatigue  contribution  for  a given  cycle,  the  crack  growth  is 
classically  divided  into  two  parts,  the  first  one  concerning  the  fatigue  propagation,  and  the 
second  one  the  creep: 


For  the  treatment  of  a general  complex  loading,  for  instance  with  hold  time  and  overload 
as  described  in  figure  1 , we  post  the  fatigue  crack  growth  (the  minimum  value  of  the  stress 
intensity  factor  is  memorised)  for  the  given  loading  at  the  next  unloading.  Only  time 
dependent  phenomena  such  as  creep  crack  growth,  threshold  relaxation  and  oxidised  zone 
increasing  are  computed  during  the  loading.  Let  us  precise  that  the  threshold  is  supposed  to 
be  constant  during  the  creep  crack  growth.  It  is  only  actualised  just  before  to  post  the 
fatigue  crack  growth,  and  just  after  for  instance  to  take  into  account  an  overload  (or  an 
under-load).  The  environment  effects  are  also  posted  for  the  next  cycle  through  decreasing 
of  the  material  toughness.  ^ 
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-4  threshold  relaxation  Kf  = f(t,  T) 
— > creep  (with  constant  threshold  Kcs) 
-4  oxidation  / = f(tt  T) 


J— > fatigue 

— » Environment ; toughness 
calculation  for  the  next 
cycle.  Kc  = /(/,) 

— » threshold  actualization 
due  to  fatigue  loading 


time 


FIGURE  1.  General  scheme  of  the  creep-fatigue-oxidation  interaction  model 


2.2  Non  isothermal  writing 

Let  us  consider  a thermo-mechanical  loading  for  which  the  temperature  T(t)  can  strongly 
vary  during  the  force  cycle  F(t)  or  the  stress  intensity  factor  K(t)  at  the  instant  t . The 
difficulty  for  the  model  extension  in  non  isothermal  condition  concerns  the  fatigue  law 
which  predicts  a fatigue  damage  rate  per  cycle.  Which  temperature  do  we  have  to  consider 
? We  can  consider  the  maximum  temperature  over  the  cycle  to  be  in  the  security  way,  or  to 
define  a mean  temperature  T'  as  proposed  by  Taira  (4). 

The  way  we  propose  here  to  account  for  the  temperature  dependency  is  to  use  a reduced 
K(t) 

stress  S = as  parameter  describing  the  thermo-mechanical  fatigue  cycle,  where 

Km(T)  is  the  material  toughness  (without  oxidation  effect)  which  is  -temperature 

dependent,  the  fatigue  law  can  then  be  supposed  temperature  independent.  In  another 
words,  this  hypothesis  requires  that  the  fatigue  material  properties  evolve  with  temperature 
only  through  the  evolution  of  the  material  toughness  supposed  here  to  obey  to  an  Arrhenius 
law: 


K„0)  = K!,ex j>(A) 


P) 


where  Q (the  activation  energy)  and  K°v  are  material  coefficients.  We  define  then  a fatigue 
law  as  a Forman  [5]  expression  of  S variable: 


c;k,-sfP 

(l-tffe-Sj 


(3) 


) Fatigue 

where  C'f  and  Tjf  are  temperature  independent  material  coefficients.  Sf  is  the  fatigue 
threshold,  also  temperature  independent.  R is  the  loading  ratio,  Sc  the  normalised 
toughness  equal  to  1 without  environment  effect,  and  SM  the  maximum  value  of  S during 
the  cycle.  From  tests  results  obtained  in  fatigue  regime,  without  time  effect  such  as  creep  or 
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environment  effect,  and  for  several  temperatures,  we  have  to  find  only  one  mistress  curve  if 


we  report 


da 

In 


as  a function  of  AS  = ■ 


K 


At 


K! 


Klv(T(t))  Kcv(T(t)) 


Tests  performed  by  Bemede  [6]  on  Astroloy  superalloy  in  vacuum  and  in  air  at  high 
temperature  show  a great  influence  of  environment  on  the  crack  propagation  process.  The 
main  environment  effect  is  the  material  oxidation  at  the  crack  tip  leading  to  a reduction  of 
the  mechanical  properties.  It  is  introduced  in  the  model  through  the  variation  of  the  local 
normalised  toughness  induced  by  oxidation: 


„ , „ K" 

( * Yl 

mz 

5c(0=^f 

cv 

1 -u  + w.exp 

Sc  = Min(S  c„  1) 


where  mz'  and  u are  material  temperature  independent  coefficients.  By  considering  the 
toughness  Kcox  ( T ) of  the  completely  embrittled  material  at  the  crack  tip,  necessarily  lower 

than  KCV(T ) and  also  supposed  to  obey  to  the  same  Arrhenius  law  (same  activation 
energy): 


(5) 

is  the  penetration 


dl  P(0  = j(F(x)a)4  l~zdt 

(6) 

a =ac  exp (^-) 

a0  and  the  activation  energy  Qo  are  material  coefficients.  This  function  traduces  the 
material  oxidation  kinetics  at  the  crack  tip  [7]  and  its  temperature  dependency.  The  function 
F(x)  accounts  for  the  loading  effect  on  the  oxidation  kinetics,  for  instance  an  hold  time 
after  an  overload,  x is  the  ratio  between  the  actual  load  and  the  previous  overload: 

F(x)=cx4=c(f^)‘  (7) 

Fm 

c is  a constant.  In  the  case  of  increasing  load,  X equals  to  1 and  the  oxidation  effect  on 
fatigue  process  is  completely  operative. 

As  creep  damage  is  concerned,  the  extension  in  non  isothermal  does  not  require 
supplementary  development  because  the  creep  model  is  described  by  a temporal  equation: 

da  = ) Cc(T(t))\K(t)  - Kl]Vmdt  (S) 

'l 

A given  cycle  (time,  load,  temperature)  is  treated  by  considering  time  steps  sufficiently 
small  to  assure  the  calculation  convergence,  and  by  taking  at  each  time  step  the 

corresponding  parameters  Cc  ( T(t ))  and  T|f  (T(t)) . Crack  propagation  by  creep  is  posted 


K ( T)  = K°  exp(-^-) 

COX  y ' COX  rVj^J»7 


we  assume  then  that  the  ratio 


K" 


is  temperature  independent.  / (t) 


K°  p 

cv 

length  of  oxide  inside  the  material  as  a function  of  time  and  temperature: 
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as  soon  as  the  actual  stress  intensity  factor  K(t)  is  higher  than  the  non  propagation 
threshold  AT/.  We  can  remind  here  that  this  threshold  is  considered  as  a constant  during  the 

cycle  and  equals  to  the  fatigue  non  propagation  threshold  K fs  calculated  at  the  previous 
cycle.  On  the  other  hand,  we  compute  during  creep  the  relaxation  of  this  one,  by  a 
relaxation  type  equation  deduced  from  the  study  of  viscous  behaviour  of  the  material: 

d(-^-)  = -A(T(t))(^-fITU)]  dt  (9) 

K(t)  K(t)J 

As  well  as  the  creep  model  integration,  we  consider  the  instantaneous  values  of  the 
parameters  A(T(t))  and  CO  (T  (t ))  at  each  time  step  of  the  calculation. 

3 APPLICATION  TO  A NICKEL-BASE  SUPERALLOY  FOR  TURBINE  DISK 

We  present  in  this  section  the  application  made  on  N18  Nickel-base  superalloy  at  high 
temperature.  A wide  experimental  program  has  been  conducted  by  SNECMA  Moteurs  [8] 
from  450°C  to  650°C  for  various  loading  conditions.  After  showing  the  results  of  the 
identification,  we  give  predictions  of  many  isothermal  and  non  isothermal  tests. 


3. 1 Identification  of  the  model  parameters 

The  proposed  formulation  makes  the  model  attractive  for  its  identification  on  a large 
temperature  domain.  As  a matter  of  fact,  the  fatigue  law  evolves  with  temperature  through 

the  variation  of  the  material  toughness  KCV(T ) supposed  to  be  known.  Fatigue  tests  results 
performed  at  high  temperature  but  at  high  loading  frequency  in  a pure  fatigue  regime,  for 
which  time  effects  such  as  creep  and  oxidation  can  be  neglected;  or  at  low  temperature 
whatever  the  loading  frequency  (no  time  effect),  are  then  used  to  identify  the  fatigue 
coefficients. 


FIGURE  2.  Identification  of  the  pure  fatigue  model 
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The  figure  2 gives  the  result  of  Cf  and  identification  (the  initial  non  propagation 

threshold  S{’  is  taken  equal  to  0,5)  by  reporting  in  a diagram  (-^-,A S)  tests  results 

dN 

obtained  at  450°C  (we  observed  any  time  effect  at  this  temperature),  550°C  (we  suppose  to 
have  a small  time  effect  for  those  tests  performed  at  0,5  Hz)  and  650°C  (at  20  Hz  in  a pure 
fatigue  regime)  for  the  same  loading  ratio.  At  a given  temperature,  we  can  then  calculate 

CAT)  and  K{  (T)  parameters  by  the  following  relations: 


Cf(T)=C)-K]<(T) 
K{  (T)  = Sf  ■ KCV(T) 


(10) 


Pure  creep  tests  at  different  temperatures  are  required  to  identify  the  creep  model 
temperature  dependent  parameters  Cc(T)  and  T)ir(7’)of  equation  8.  The  initial  non 


propagation  threshold 


is  also  taken  equal  to  0,5 . Wc  reported  in  figure  3 the  results 


of  the  identification  at  650°C. 


FIGURE  3.  Identification  of  the  pure  creep  model  at  650°C 


To  complete  the  information  obtained  by  tests,  all  the  phenomena  concerning  the  crack 
opening  and  closure  and  the  stress  redistribution  at  the  crack  tip  induced  by  plasticity  or 
visco-plasticity  during  the  crack  propagation  can  be  obtained  by  a numerical  analysis  using 
a finite  element  method,  as  it  has  been  made  by  Prigent  (1)  in  the  case  of  Astroloy 
superalloy.  This  numerical  analysis  conducts  to  the  identification  of  a(T)  and 
0(7’)  threshold  relaxation  parameters  of  equation  9 (the  visco-plastic  behaviour  of  the 
material  at  high  temperature  is  of  course  supposed  to  be  known  and  predicted  by  accurate 
constitutive  models). 


The  oxidation  kinetics  of  the  material  as  well  as  its  temperature  dependency  is  (aoand 
Q0  coefficients  of  equation  6)  supposed  to  be  known.  The  evolution  of  the  material 


toughness  ahead  of  the  crack  tip  ( K°cox  and  then  the  ratio 


K". 


; mz‘  and  u coefficients  of 
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equation  4)  can  be  deduced  from  special  tests  performed  on  oxidised  specimens.  The 
coefficient  c (equation  7)  which  traduces  the  loading  effect  on  the  oxidation  kinetics  is 
identified  by  one  test  result  obtained  in  fatigue-creep  regime,  for  instance  by  the  simulation 
of  a test  with  hold  time. 

3.2  Predictions  of  isothermal  tests 

We  give  in  this  section  several  comparisons  model/prediction  of  tests  performed  on 
specimens  with  notches  of  different  geometry  (KB,  KBR  and  FC). 


FIGURE  4.  Prediction  of  time  and  load  ratio  effects  at  550°C 

In  figure  4 are  reported  the  model  predictions  at  550°C.  As  it  can  be  seen,  the  model 
predicts  almost  any  time  effect  at  this  temperature  in  agreement  with  the  tests  results.  The 
load  ratio  effect  is  also  correctly  reproduced.  Complex  loading  have  been  defined  to 
simulate  on  a specimen  what  can  be  observed  on  the  real  component  such  as  SF1,  SF2  and 
A’20  cycles  reported  on  figure  5. 


FIGURE  5.  Complex  cycles  SF1,  SF2  and  A’20 

The  figure  6 gives  the  model  predictions  for  these  complex  cycles.  The  loading  rate  effect 
is  well  reproduced  as  SF1  and  SF2  cycles  are  concerned.  The  test  reproducing  A’20 
mission  is  predicted  with  a little  too  high  crack  propagation  rate. 
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FIGURE  6.  Predictions  of  complex  tests  SF1,  SF2  and  A’20  at  550°C 


We  present  now  tests  predictions  at  650°C  for  which  time  effects  are  important. 


FIGURE  7.  Predictions  of  tests  with  hold  time  at  650°C 

In  figure  7 are  given  the  predicted  curves  for  tests  with  different  hold  times.  The  model 
predictions  are  quite  good  whatever  the  hold  time  duration.  It  appears  a small  difference 
between  the  two  predicted  tests  with  hold  times  of  300  s and  900  s,  in  agreement  with  the 
tests  results. 
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FIGURE  8.  Predictions  of  tests  with  hold  time  for  different  load  ratios  at  650°C 

The  last  predictions  reported  on  figure  8 concern  the  loading  ratio  effects  at  650°C  for  tests 
with  the  same  hold  time  of  300  s.  The  loading  ratio  effects  predicted  by  this  crack 
propagation  model  reproduce  quite  correctly  what  it  can  be  observed  with  N18  superalloy 
and  for  tests  with  consequent  creep-fatigue-oxidation  interaction. 

3.3  Predictions  of  complex  tests  under  non  isothermal  condition 

Three  different  complex  loading  in  non  isothermal  condition  were  studied  [8]  to  simulate 
on  specimens  what  can  be  subjected  to  critical  zones  of  a turbine  disc. 


FIGURE  9.  Thermo-mechanical  fatigue  cycles  simulated  in  the  laboratory 
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FIGURE  1 0.  Comparisons  between  calculated  and  experimental  lifetimes 

The  applied  load  to  the  specimen  and  the  temperature  loading  are  reported  in  figure  9. 
The  temperature  varies  for  each  cycle  from  about  200°C  to  600°C.  The  lifetime  is  the 
number  of  missions  to  reach  a given  crack  length.  We  present  in  figure  10  the  comparisons 
between  the  experimental  and  calculated  lifetimes  for  isothermal  and  non  isothermal  tests. 
The  model  predictions  are  fairly  good,  as  well  as  for  many  isothermal  tests  at  several 
temperatures  and  loading  conditions,  and  for  the  thermo-mechanical  fatigue  tests  for  which 
the  lifetime  is  correctly  predicted,  even  if  the  model  parameters  have  not  been  identified  at 
the  maximum  temperature  of  the  cycles  (600°C)  but  interpolated  with  results  obtained  at 
550°C  and  650°C. 

4 CONCLUSION 

The  works  summarised  in  this  paper  demonstrate  actual  possibilities  to  predict  the  crack 
propagation  of  metallic  materials  at  high  temperature  such  as  Ni-base  superalloy  for  turbine 
discs.  The  strong  requirements  in  design  procedures  have  led  to  the  development  of 
sophisticated  crack  propagation  models  to  account  for  fatigue,  creep  and  environment 
effects  as  well  as  theirs  interactions  on  crack  propagation  rate.  We  proposed  a 
phenomenological  fatigue-creep-oxidation  interaction  model  for  an  application  in 
isothermal  condition  but  also  in  non  isothermal  condition.  We  presented  the  lifetime 
prediction  of  many  isothermal  tests  and  of  a few  anisothermal  tests.  Predictions  are  fairly 
good  for  the  set  of  experimental  data  and  are  very  promising  to  make  the  model  predictive 
under  complex  thermo-mechanical  loading. 
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